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Abstract 

The leaky box model is now ruled out by measurements of a cosmic ray gradient 
throughout the galactic disk. It needs to be replaced by a more refined treatment 
which takes into account the diffusion of cosmic rays in the magnetic fields of the 
Galaxy. We have estimated the flux of antiprotons on the Earth in the framework of 
a two-zone diffusion model. Those species are created by the spallation reactions of 
high-energy nuclei with the interstellar gas. Another potential source of antiprotons 
is the annihilation of supersymmetric particles in the dark halo that surrounds our 
Galaxy. In this letter, we investigate both processes. Special emphasis is given to 
the antiproton signature of supersymmetric dark matter. The corresponding signal 
exceeds the conventional spallation flux below 300 MeV, a domain that will be 
thoroughly explored by the Antimatter Spectrometer experiment. The propagation 
of the antiprotons produced in the remote regions of the halo back to the Earth 
plays a crucial role. Depending on the energy, the leaky box estimates are wrong 
by a factor varying from 0.5 up to 3. 
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1 Introduction 



The flux of antiprotons on the Earth is about to be measured with unprecedented accuracy 
by the Antimatter Spectrometer experiment (AMS) [p]]. A detector will be shipped in 
orbit on the space shuttle. Cosmic rays will be directly detected in space. Unlike in 
balloon-borne experiments, measurements will not be spoilt by the atmosphere. Even in 
the stratosphere, the column density of air is still of the order of a few g-cm -2 . This value 
is still an order of magnitude above the average of the column density of interstellar gas 
toward the region of the galactic centre. The AMS mission will be mostly sensitive to 
the energy range extending between 100 MeV and a few GeV, in a domain so far barely 
explored and where previous measurements are inconsistent with one another. The proton 
and antiproton fluxes will give precious indications on the propagation of cosmic rays at 
low energy. 

As detailed in section 2, antiprotons are predominantly produced by the interactions 
of primary cosmic rays, mostly protons, with the interstellar gas that pervades the thin 
matter disk of our Galaxy. Quite exciting is the possible presence of an additional an- 
tiproton flux at low energy, which would signal the presence of massive neutral particles 
around our Galaxy. The Milky Way is embedded inside a halo of invisible material whose 
nature has not yet been determined. That dark matter could be made of relics from the 
big-bang. In that case, supersymmetric species are a favoured option and AMS should be 
sensitive to the antiprotons produced by their mutual annihilations. 

Our understanding of cosmic ray propagation throughout the Galaxy has recently 
improved. In the old leaky box scheme, the Galaxy is pictured as one single region where 
production, interactions and diffusion of cosmic rays take place homogeneously. Cosmic 
rays are stored in a vast domain, from which they escape after being confined for some 
time. According to the leaky box model, they are uniformly distributed. This last feature 
is not quite correct. Measurements of the 7-ray diffuse emission of the Galaxy actually 
point toward a gradient of the cosmic ray density along the galactic ridge ||. The proton 
flux varies with a galactocentric radius r. Recently, the 7-ray diffuse radiation of the 
Perseus arm region has been accurately determined by the EGRET experiment on board 
the Compton Gamma Ray Observatory. The 7-ray emissivity of interstellar hydrogen 
atoms immersed in the cosmic radiation drops by a factor of ~ 1.7 between the solar 
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circle (r ~ 9 kpc) and the Perseus arm (r ~ 13 kpc) [Q. The naive leaky box model needs 
actually to be replaced by a more refined treatment of cosmic ray propagation. 

In section 3, a two- zone diffusion model is presented. Cosmic ray primaries are pro- 
duced inside a thin matter disk in which they propagate and interact with the interstellar 
gas. Those spallations generate antiprotons, among a variety of secondary nuclei. Cosmic 
rays diffuse in the chaotic magnetic fields carried by the galactic stellar winds. That dif- 
fusion takes place in an extended region, which encapsulates the previous thin disk and 
extends a few kiloparsecs away from the galactic plane. In this letter, we focus on the 
antiproton flux on the Earth, which is due to the annihilation of supersymmetric particles 
in the galactic halo. The corresponding antiproton yield is also derived in that section, 
in the framework of diffusion. 

Finally, in section 4, we discuss our results. Massive neutral species annihilate in 
remote regions of the Galaxy, and antiprotons propagate back to the solar system. Their 
diffusion from the outer parts of the halo all the way toward the Earth is a crucial 
mechanism, which cannot be accounted for by the naive leaky box model. However, 
the main trend is preserved: the supersymmetric antiprotons show up at low energy. 
Above ~ 1 GeV, they are swamped in the flux of spallation antiprotons. The strongest 
supersymmetric signal is generated either by a pure gaugino or a pure higgsino. The 
intermediate case, where the supersymmetric species is a mixture of a gaugino and a 
higgsino, will be barely detectable, unless the abundance of that particle in the galactic 
halo is larger than its cosmological average. We also compare our predictions to the leaky 
box estimates. A correct treatment of the cosmic ray diffusion throughout the Galaxy is 
important. In the leaky box model, the antiproton flux is underestimated at low energy 
by a factor of 2 and overestimated at high energy by a factor of ~ 3. We have also derived 
the confinement time of the antiprotons inside the galactic halo and disk as a function of 
energy. Such a result may be directly incorporated in a leaky box model to yield a not 
really bad estimate of the antiproton flux. We feel that it may be helpful to colleagues 
who would like to rapidly derive a fair value of the antiproton signal without getting into 
the intricacies of our more involved treatment. 
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2 The sources of antiprotons 



Antiprotons are naturally produced by the nuclear interactions of cosmic ray particles 
with the interstellar medium. In order to estimate the energy spectrum of the antipro- 
tons created during a single collision, we have used the Lund Monte Carlo simulations 
PHYTIA and JETSET [|]. Since both cosmic rays and the interstellar medium are mostly 
composed of protons, the fundamental process is the collision of a high-energy proton with 
a proton at rest in the interstellar medium. This reaction produces hadrons, mostly pions 
whose radiative decays are the principal source of the galactic diffuse gamma-ray emis- 
sion. During the fragmentation process, when the inelastic collision produces coloured 
strings, resonances such as A , A~, A or S~ are formed and subsequently decay into 
antiprotons. Spallation reactions also produce antineutrons, which in turn decay into 
antiprotons. The threshold for the antiproton production mechanism p + p — > 3p + p 
is 7 GeV. The energy E p of the incoming cosmic ray proton has been varied from that 
threshold up to 2.7 TeV. For each value of E p , a million events have been generated to 
yield the antiproton multiplicity per collision 

Y p (E p -> E p ) = ^ , (1) 

in the antiproton momentum range extending from 100 MeV to 100 GeV. 

One of the best motivated particle candidates for cold dark matter is provided by 
supersymmetry ||. The lightest supersymmetric particle (LSP) is stable, provided R- 
parity is conserved. In a large region of the supersymmetric parameter space, that species 
is the lightest neutralino (x), defined as the lowest-mass superposition of the gaugino 
(7, Z) and the higgsino (Hi,H 2 ) fields : 

X = Oi7 + a 2 Z + a 3 Hi + a 4 H 2 . (2) 

The possibility to detect the presence of neutralinos in our galactic halo has been studied 
extensively in many different ways ||. Direct detection of the nuclear recoil induced 
by neutralino-nucleus elastic scattering is, both theoretically [J7|] and experimentally ||, 
under investigation. The indirect detection of high-energy neutrinos coming from the 
centre of the Earth or the Sun ||, or the indirect detection of exotic components in the 
primary cosmic rays, such as 7, e + and p, have been studied in detail |TI], |nj. Here we 
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examine the latter possibility, and concentrate on the diffusion of the cosmic antiprotons 
p from the remote regions of their production back to the inner Galaxy and the Earth. 

The differential rate (per unit volume and per second) for the production of antiprotons 
from x~X annihilations is given by 



-y(E p ) = (av)f(E P ) \*L , ( 3 ) 



m x 



in which Ep denotes the antiproton energy, a is the x~X annihilation cross section and v 
is the neutralino velocity in the galactic halo. The neutralino density p x varies with its 
position. For a single annihilation, the antiproton energy spectrum is 

/W) . ( I } {*te^±S}. s ^ @, (4) 



a 



F,f \ 



where F describes the x~X annihilation final state and B^J is the branching ratio into 
the quarks or gluons / in the channel F. The differential distribution of the antiprotons 
generated by the hadronization of quarks (with the exception of the top quark) and of 
gluons is denoted by dNp / dEp and depends on the nature of the species /. Of the various 
quantities that are present in Eq.(^), (av) and f(Ep) depend on the neutralino properties. 
The antiproton production rate also depends on the distribution p x of the dark matter 
particles inside the galactic halo. At this stage, some comments are in order : 

(i) The neutralino annihilation cross section (av) is evaluated as in Ref. Jy], with the 



parameters of the minimal supersymmetric standard model (MSSM) fixed at the following 
values: the lightest Higgs mass is = 55 GeV while the sfermion mass is fh = 500 GeV 
and tan/3 = 8. The x composition parameter P = a\ + a\ is fixed, for the three different 
cases considered here, at the values P = 0.01 (higgsino), P = 0.5 (mixture) and P = 0.99 
(gaugino) . 

(ii) For the p differential distribution f(Ep), we have evaluated the branching ratios Bj 
for all annihilation final states that may produce antiprotons, i.e. direct production of 
quarks and gluons, generation of quarks through the intermediate production of Higgs 
bosons, gauge bosons and the top quark. The distributions dNp / dEp from the hadroniza- 
tion of quarks (with the exception of the top quark) and gluons have been computed by 
using the Monte Carlo code JETSET 7.2 11. 
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(iii) The neutralino halo distribution p x is taken to be spherically symmetric. In the 
axisymmetric coordinate system r and z, the density profile is given by 



p x (r,z) = p x (Q) j-A±ZJ_j j (5) 



where a = 3.5 kpc is the core radius of the dark matter halo. Particular care must be 
taken about the local neutralino density p x {&), which depends on the LSP properties. 



Actually, if the big-bang relic density Q x h 2 , which we evaluate following fl2|], is too small 



to account for the cosmological dark matter, the density of neutralinos in the galactic halo 
should be corrected by a factor of r] x . The latter deals with the fact that the neutralino 
density is less than the halo density whenever Q x h 2 is smaller than a minimal value of, 
say, (Qh 2 ) min = 0.03, which is compatible with the observed rotation curve of the Galaxy. 
At cosmological distances, this ratio is given by 



n x h 2 



/ n in i 



whereas in the galactic halo, it may well be significantly larger. For instance, there may 
be segregation between neutralinos and dark baryons as a result of the dissipation of the 
latter. Therefore, we must keep it in mind that r] x could be as small as its cosmological 
average, but it may also be way larger. 



3 The diffusion model for antiprotons 

Parker has studied the propagation of cosmic rays inside the Galaxy as a consequence 
of their scattering by the irregularities of magnetic fields. The presence of the latter is 
now firmly established by synchrotron radiation far above the plane of our Galaxy as 
mentioned by Badwar Magnetic fields are also detected in other galaxies [pi}] . So 
the cosmic ray transport in the Galaxy crucially depends on the diffusion across magnetic 
fields. In the following, we will assume an isotropic diffusion with an empirical value for 
the diffusion coefficient. Another proof of the existence of a diffusion region in our Galaxy 
is provided by the recent diffuse 7-ray observations, which point toward the presence of 
cosmic rays far above the galactic disk. As regards convection, it seems probable that 
cosmic rays as well as the stellar wind from disk stars contribute to push out the magnetic 
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fields so that the region where the particles diffuse inflates at a speed of the order of a 



few kilometers per second [13, 15]. Because convection has been shown to be negligible 
TEfl , we will disregard it and will focus our analysis on the pure diffusion case. Thus, 



our Galaxy can be reasonably modelled by a thin disk of atomic and molecular hydrogen 
with < r < R = 20 kpc and \z\ < h = 100 pc, associated to an extended region of 
diffusion containing irregular magnetic fields with the same radial extension and \z\ < L 
= 3 kpc. Those various regions are superimposed to the spheroidal halo of dark matter 
whose density profile has already been given in relation (|5|). That two-zone diffusion model 
is in good agreement with the observed primary and secondary nuclei abundances Jl( 



Since their discovery by Golden in 1979 [17], cosmic antiprotons have been thoroughly 
studied in the framework of the leaky box model [18|. We analyse here their propagation 
throughout the Galaxy in the light of a two-zone diffusion model. 

In a stationary regime, the propagation equation of cosmic antiprotons may be ex- 
pressed as 

Qti- 

— = = KAn p - T p n p + q p , (7) 

where np(Ep,r, z) is the density of antiprotons with energy E p at location (r, z). In the 
right-hand side of relation (0), the first term describes the diffusion of antiprotons. The 
diffusion coefficient K is constant at low energies, but raises with the rigidity p of the 
particle beyond a rigidity of 3 GV. This behaviour can be modelled by the form 

K = 6 x 10 23 m 2 ■ s- 1 (l + • (8) 

The second term in Eq.(]7J) describes the destruction of antiprotons by their interactions 
with the interstellar medium. That destruction rate is shown to be very small but has 
not been neglected in what follows. The collision rate of antiprotons with the interstellar 
hydrogen is given by 

Tp = Opu v p n H , (9) 

where <jpH is the total antiproton interaction cross section with protons, Vp denotes the 
velocity and nn = 1 cm -3 is the average hydrogen density in the thin matter disk. The 
last term in relation (|7J deals with the various sources of antiprotons. Those species are 
produced by the spallation of cosmic protons on the interstellar matter of the disk. The 
antiproton production rate involves a convolution over the incident cosmic proton energy 
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spectrum dn p /dE p 

qf{E p ) = J dE p |^r(£p)j T p (E p )Y p (E p ^E p ) . (10) 

The collision rate T p of protons with the interstellar gas is defined in just the same way as 
the collision rate T p of antiprotons in relation @. The antiproton differential spectrum 
Y p produced during a single proton spallation was previously defined in Eq.([l|). The 
super symmetric source term (|3]) may also be present if antiprotons are produced by the 
mutual annihilation of neutralinos. This leads to an additional antiproton flux. 

We have solved Eq.(^) following the analysis by Webber, Lee and Gupta RIB] . At 
the edge of the domain where cosmic rays are confined, the particles escape freely and 
the diffusion becomes inefficient. Thus the density vanishes at the boundaries of the 
domain where cosmic rays are confined by diffusion. This provides the initial conditions 
for solving the diffusion equation. Then, because the problem has a cylindrical symmetry, 
the densities n p and n p may be expanded as series of the Bessel functions of zeroth order 
Jo (Ci x ) where Q is the ith zero of Jo and where x = r/R. Details may be found in fl6| . 
The cosmic ray sources are located in the thin gaseous disk and we have taken the radial 
distribution of supernovae remnants and pulsars measured by Lyne et al. [fil| . As regards 
the spallation mechanism, suffice it to say that an effective antiproton multiplicity may 
be defined as 

Yf(E f ) = j;~ d E p {l^}Y f{ E p ^E f ) , (11) 

so that no convolution with the cosmic proton energy spectrum is needed any longer. The 
differential flux of protons of energy E p is denoted here by $ p (E p ). 

In the case of the supersymmetric antiprotons, the resolution of Eq.(0) is more in- 
volved. The supersymmetric source term q p nsy now depends on both the galactocentric 
radius r and on the vertical coordinate z. In the solar neighbourhood, the antiproton 
energy spectrum due to the mutual annihilations of exotic particles is found to be : 

dn p ^ f 4 1 f T ( ^ R Q 



dE p 



;(©) = E {j-} { J o(o±f)} K 2 ^} {/o ^)QT y (z)dz\ , (12) 
where the Bessel transform Q^ usy (z) is defined as 

QT Sy (z) = £ x dx q^(E p , r = xR, z)J (Ct^jf) ■ (13) 
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The vertical distribution Ti{z) is given by 

sin! 

= ^\c. ' , (14) 



sinh | 


l 2 




.)} 


sinh | 


[14 





where the parameter Si is equal to 2Q/R. Finally, the coefficient Ai in relation (|T|) stands 
for 

= 2hT p + ATSi coth j ^lJ , (15) 

and only depends on the antiproton energy Ep through the diffusion coefficient K (|§D and 
the collision rate Tp @. Expression ([12]) involves an integral on both the galactocentric 
radius r and on the height z. That integration has been performed here with the specific 
form (|3|) for the supersymmetric production rate qp usy . 



4 Results and discussion 

We first show, in fig. 1, the standard interstellar p/p ratio obtained in the framework of our 
diffusion model, as a function of energy. As there is some uncertainty in the measurement 
of the spectral index a of the high-energy proton spectrum, we considered the cases 
a = 2.65 (dashed line), 2.70 (intermediate solid curve) and 2.75 (dotted line). The p/p 
ratio decreases at low energy because the low-energy antiprotons must be produced with a 
large backward momentum in the centre-of-mass reference frame, and so their progenitors 
are very high energy protons, whose density is very low. The harder the spectrum, the 
larger the p/p ratio. Our results are in fair agreement with the estimates of Gaisser and 
Schaefer |Tj| collected in their fig. 5. Differences arise from the improved treatment of 
the diffusion mechanism with respect to the phenomenological leaky box model. 

The antiproton signature p/p of neutralino pair-annihilations in the halo must be 
compared with the standard p/p ratio discussed above. Figure 2 features this ratio for 
three species of neutralinos, as discussed in section 2. At low energy, the SUSY antiproton 
production may exceed the standard signal in the pure gaugino and higgsino cases. In 
the intermediate situation, it is well below the secondary antiproton flux because the 
neutralino relic density is Q x h 2 ~3x 10~ 3 , only a tenth of the minimal value (Qh 2 ) m i n . 



8 



This is not the case for a pure gaugino or higgsino. Note that the suppression factor r\ x 
may be larger than 1/10 so that the dotted curve may well be shifted upward. 

In the leaky box model, cosmic rays are distributed homogeneously over the entire 
diffusion domain where they are confined by the magnetic fields of the Galaxy. Their 
diffusion may be described from a phenomenological point of view by the confinement 
time r, which measures the rapidity with which the particles manage to escape in the 
intergalactic medium. That confinement duration r is related to the column density, i.e. 
the amount of matter which the cosmic rays cross during their erratic journey. Typical 
values of a few g-cm~ 2 are necessary to convert primary nuclei such as carbon, oxygen 
or nitrogen into secondaries such as boron or beryllium. This implies a confinement time 
in the matter disk of the Galaxy of order 10 7 years. Because the unstable isotope 10 Be 
decays with a half-life of 1.6 megayears (My), it provides a unique chronometer of the 
time actually spent by cosmic rays since their production. Measurements of its abundance 
relative to its stable partner indicate that beryllium is confined during ~ 10 8 years, hence 
the need of an extended region of diffusion which encapsulates the thin gaseous disk. 
In the leaky box model, the supersymmetric production rate g| usy is averaged over the 
confinement region and multiplied by the time r g it takes for a high-energy particle to 
escape from that diffusion domain 



M/(30 >r g . (16) 



The square of the neutralino density is easily averaged over the diffusion box. More un- 
certain is the determination of the confinement time of antiprotons. Previous calculations 



fTTJ, [L1J] are based on a typical escape time r g of order 10 years. The authors nevertheless 
note that estimates for that value may well vary by an order of magnitude. Our two- zone 
diffusion model directly incorporates that escape time, which depends on the cosmic ray 
distribution throughout the Galaxy and on the energy through the diffusion coefficient. 
The escape time from the disk and from the whole Galaxy r g may also be derived 
from that model. For protons or antiprotons, the net production rate Q in the entire 
Galaxy may be computed. A fraction of the particles that are produced are destroyed 
by spallation on the interstellar gas of the disk. What remains is the escape rate N esc , 
provided that steady state is achieved 

Q = TN d + N csc . (17) 
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The confinement time of cosmic rays in the thin gaseous disk and in the whole domain of 
diffusion are related to the escape rate through 



"1 = tli. 



N = — = 

1 v esc 



where the total number of particles in these two regions are respectively denoted by 
and N g . For both protons and antiprotons, we have evaluated the various escape times 
in the framework of our diffusion model. We find that they are fairly similar. The solid 
curve in fig. 3 features the variations of the antiproton escape time Tg clS cL function of 
kinetic energy. The dashed line stands for the escape time from the matter disk. Both 
escape times decrease with energy. As a matter of fact, the diffusion coefficient increases 
at high rigidity so that particles diffuse, and therefore escape, more easily. From 100 MeV 
to 1 GeV, the disk confinement duration is ~ 11 My while the escape time from the 
Galaxy as a whole is r g ~ 180 My, a factor of 2 above the value currently used in the 
previous estimates of the antiproton flux. In fig. 4, we focus on the antiproton signal 
of supersymmetric dark matter. For purposes of comparison, the ratio of the leaky box 
prediction fll6|) to the exact result (0) is plotted as a function of the cosmic ray kinetic 
energy. The solid line corresponds to an escape time r g of 3 x 10 15 s, independent of 
the energy. At low energy, the leaky box model underestimates the antiproton flux by a 
factor of 2. On the contrary, it is too optimistic at high energy by a factor ~ 3. If in 
relation (|l6|) , we use the value of the escape time T g , which has been derived in fig. 3, the 
situation considerably improves as featured by the dashed curve. The leaky box model 
may actually be used, provided the escape times are correctly derived with a diffusion 
model. The excess over relation ( |T2"D is at most 30% for a kinetic energy of 1 GeV. 

In this letter, we have used a two-zone diffusion model to describe the propagation of 
antiprotons throughout the Galaxy. We have estimated the conventional and supersym- 
metric antiproton yields on the Earth. As regards the neutralino pair-annihilations in the 
halo, a correct treatment of diffusion improves the old leaky box estimates. We have pre- 
sented here three characteristic examples of neutralino composition. A more exhaustive 
investigation of the supersymmetric parameter space is now necessary. 
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Figure Captions 



Fig 1 : The p/p ratio obtained in the framework of a diffusion model is displayed as a 
function of energy. The sources of antiprotons are the spallation reactions of cosmic ray 
primaries on the interstellar gas. The spectral index of the proton flux has been varied 
from 2.65 (top curve) to 2.75 (bottom curve). 

Fig 2 : The p/p ratio is featured as a function of energy, for various antiproton sources: 
standard spallations (solid line) assuming a spectral index of 2.7 for the primary proton 
distribution, and neutralino pair annihilations, for three different neutralino compositions 
(gaugino, higgsino, and a mixture). 

Fig 3 : The diffusion time scales of antiprotons in the thin gaseous disk (dashed line) 
and in the whole Galaxy (solid line) are plotted as a function of energy. Those time scales 
are inferred in the framework of the diffusion model discussed in the text. 

Fig 4 : The ratio of leaky box results over direct diffusion estimates are presented as 
a function of energy, for the antiprotons from neutralino pair annihilations. The solid 
line is obtained by setting in the leaky box model the antiproton escape time scale at the 
constant value of 3 x 10 15 s. In the case featured by the dashed curve, this escape time 
scale has been borrowed from fig. 3. 
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Figure 3 
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